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ABSTRACT 

We compare standard models of accretion disks around black holes that include the ap- 
propriate zero-torque inner boundary condition and relativistic effects on the emission 
and propagation of radiation. The comparison is performed adopting the multicolor 
disk blackbody model (MCD) as reference and looking for the parameter space in 
which it is in statistical agreement with "more physical" accretion disk models. We 
find simple 'recipes' that can be used for adjusting the estimates of the physical in- 
ner radius of the disk, the black hole mass and the accretion rate inferred using the 
parameters of the MCD fits. We applied these results to four ULXs for which MCD 
spectral fits of their X-ray soft spectral components have been published and find that, 
in three cases (NGC 1313 X-1, X-2 and M 81 X-9), the black hole masses inferred for 
a standard disk around a Schwarzschild black hole are in the interval ^ 100 — 2OOAf0. 
Only if the black hole is maximally rotating are the masses comparable to the much 
larger values previously derived in the literature. 
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1 INTRODUCTION 

Studying the spectra of X-ray binaries is one of the main 
tools for understanding the physics of these source and de- 
termine the properties of the acccretion disk and the com- 
pact object that they host. From an observational stand- 
point, Galatic black hole candidates (BHCs) show a number 
of spectral states, defined in terms of the components in 
their X-ray spectra and the flux level (e. g. [McClintock fc| 
[Remillard 2056] [Remillard fc McClin tock"2006[) . In t he so 
called high state ( [Tananbaum et al. |1972, ,Oda 19771, the 
spectrum is dominated by a soft, thermal component up to 
1-2 keV, with a power-law emerging at higher energies. In 
the low state ( Oda et al.|1971 Tananbaum et al.|1972 I the 
soft component is not present and the spectrum is well fitted 
by a power-law with a cut-off at few tens of keV. Transitions 
to an intermediate state (Mendez & van der Khs|1997| [BeT 
|loni et al.|1997p or a very high state ( [Miyamoto et al.|1991[ ) 
are also observed in which both components are present and 
equally important in terms of the emitted fiux. It is worth 
mentioning that this classification is based also on the tim- 
ing properties of the sources (e.g. Mendez et al. 19981. It 



is largely accepted that the soft, thermal component origi- 
nates from an accretion disk, while the power-law is likely 
produced by Comptonization of the soft disk photons from 
a hotter phase of the accretion flow (usually referred to as 
corona). In these assumptions, in principle the properties of 
the disk and the physical parameters of the accreting black 



hole (BH) can be inferred from the study of the soft com- 
ponent. 



Similar spectral states (e.g. Foschini et al. 2002 Feng 



fc Kaaret|2005 1 and correlated spectral/fiux variability (e.g 



La Parola 2001 Mucciarelli et al. 20071 are also observed 
in Ultraluminous X-ray Sources (ULXs), very bright, point- 
like X-ray sources in the field of nearby galaxies (see, e.g.. 



Fabbiano||2006 \. In fact, several pieces of evidence indicate 
that the majority of ULXs are X-ray binaries, chiefly among 
them the detection of periodic modulations in the X-ray flux 
(e.g. Kaaret et al.|20()6 Fabbiano et al.|2006 1 and the iden 



2004 Kaaret et al. 2004 Mucciarelh et al. 2005 1 . The spectra 



tification of stellar optical counterparts (e.g. Liu et al.|2002 



of several ULXs show a soft, thermal component, similar to 
that observed in BHCs, which has been interpreted as emis- 
sion from an accretion disk, although for some ULXs with 
spectral curvature above 2-3 keV different spectral models 
have also been proposed (e.g. Stobbart et al.|2006 Gonalves 



fc Soria|200"6} [Mizuno et al.|2007 l 

To model the soft component in BHCs and ULXs the 
so called multicolor disk blackbody model (MCD; [Mitsuda] 
et al.|1984 1 has been often adopted because it is easy to use 
and is efficiently implementated in the X-ray spectral fitting 
package XSPEC. The fitting parameters of the MCD model 
(the disk inner temperature and the normalization of the 
spectrum) depend on the BH mass and accretion rate, and 
can in principle be used to estimate the physical parame- 
ters of the accretion fiow and of the BH. However, it is well 
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known that this model represents an approximation of the 
standard, Shakura & Syunyaev (1973) disk and differs from 
it for the expression of the temperature profile, in which the 
effects of the viscous torque at the inner boundary are ne- 
glected. Also, the locally emitted spectrum is assumed to be 
a blackbody, neglecting the vertical temperature structure 
and atmospheric radiative transfer effects. Although the ef- 
fects of a zero-torque inner boundary in a pseudo-Newtonian 



potential (Paczynsky & Wiita 19801 were implemented in 
the diskpn model ( Gierlinski et al.||1999 l, temperature pro- 
file correction factors for the MOD model were also derived 
(e.g. [Kubota et~aL]|1998| [Makishima et al.|[2000|). At the 



same time, limited deviations from LTE caused by radiative 
transfer effects were accounted for using an hardening factor 
/ ~ Tcoi/Teff, where Tcoi and Te// are the color and effec- 
tive temperatures (e.g. Shimura fc Takahara|19% Zampieri 



et al.|200T I. Even adopting all these correction factors, devi 



ations of the MCD with respect to a standard disk spectrum 
may be significant if the physical state of the disk changes or 
the coupling bewteen the disk and the corona becomes im- 
portant, which may occur especially at high accretion rates 
I Merloni et al.|2000| ). 

Despite all these uncertainties, the parameters of the 
MCD fit have been used for estimating the mass of the BH 
and the accretion rate in X-ray binary systems. This has re- 
cently been done also for ULXs for which, at present. X-ray 
spectroscopy provides one of the few available methods to 
obtain zero-th order estimates of the BH mass. Indeed, if 
the soft component observed in ULX spectra is fitted with a 



MCD model (e.g. |Miller et al.|200 3, 2004a), its characteristic 
temperature is typically ~ 200 eV, much lower than that ob- 
served in Galactic BHCs. Both the low temperature and high 
normalization constant of the MCD model have been used to 
estimate the BH masses, typically obtaining values largely 
in excess of 100 M© . Modelling the soft component in these 
terms is clearly possible if ULXs are in an accretion stage 
similar to that of BHCs in the soft (or nearly soft) state, 
so that a disk spectral component is present. Along with 
the very high (isotropic) luminosity of these sources, these 
"spectroscopic estimates" have led to the suggestion that 
they may contain intermediate mass black holes (IMBHs). 

In this paper we revisit and compare standard mod- 
els of accretion disks around black holes that include the 
appropriate zero-torque inner boundary condition and/or 
relativistic effects on the emission and propagation of ra- 
diation. They are derived in the usual assumptions that: 
(1) the accretion disk is optically thick, geometrically thin 
and in a quasi-steady state; (2) the locally emitted spectrum 
from a small annulus of the disk is a blackbody {Iv = Bv), 
with a color correction factor / taking into account for atmo- 
spheric transfer effects. The comparison is performed adopt- 
ing the MCD as reference model and looking for the parame- 
ter space in which it turns out to be in agreement with "more 
physical" accretion disk models. The fits are performed keep- 
ing either the BH mass or the accretion rate fixed and al- 
lowing the inner radius of the MCD to vary. This gives us 
simple 'recipes' that can be used for adjusting the estimates 
of the physical inner radius of the disk, the black hole mass 
and the accretion rate inferred using the parameters of the 
MCD fits. 

The paper is organized as follows. Section 2 contains a 
short summary of the basic structure equations of the accre- 



tion disk models considered in this work; Section 3 presents 
the results of the comparison among the different spectral 
models; Section 4 is devoted to an application of the re- 
sults of this comparison to some representative ULX spectra, 
while Section 5 summarizes our conclusions. 



2 ACCRETION DISK SPECTRA 

In this section we summarize the basic equations of the ac- 
cretion disk models that we have considered in our work. We 
start with recalling the main properties of the standard disk 
model I Shakura fc SyunyaevHl973) and the m ulticolor disk 
blackbody model (MCD, Mitsuda et al.|1984 l, that are usu- 
ally adopted in fitting the spectra of X-ray binaries. Then, 
we shortly review the structure equations for a relativistic 
accretion disk orbiting around a Schwarzschild and a Kerr 



black hole (Bardeen et al. 1972 Novikov & Thorne 1973 



Page fc Thorne|1974| ). The emergent spectrum is calculated 
taking into account the relativistic effects on the emission 
and propagation of light (gravitational redshift, light bend- 
ing, Doppler shift) using the formalism of the transfer func- 



tion introduced by Cunningham (19751. More sophisticated 



treatments of photon propagation in curved space-times. 



both numerical (Laor|1991';'Karas et al. 


1992 Bromley etal. 


1997 


Fanton et al.||1997( Wilms„1998,, 


Ranch fc Blandford 


1994 


MuUer fc Camenzind 


2004| Beckwith fc Done |2004 ) 


and also analytical (Cadez fc Gomboc|1996 Gomboc||2001 


Cadez et al. 2003 Calvani 


20041 are available. Such tech- 



niques adopt ray tracing of the photon geodetics from the 
disk plane to the observer and are particularly suited when 
dealing with the calculation of spectral line profiles. How- 
ever, we are essentially interested in comparing the contin- 
uum spectral shape of the MCD disk with that of a standard 
and relativistic disk and, to this end, a treatment of the rel- 
ativistic geometrical optics a la |Cunningham| ( |1975[ ) appears 
acceptable. 



2.1 Standard accretion disk and MCD 



Let's consider the Standard accretion disk model (Shakura 
fc Syunyaev|1973 1 in the usual approximations (Newtonian, 
optically thick, geometrically thin and in a quasi-steady 
state). The expression of the temperature profile is obtained 
setting the energy rate of viscous dissipation equal to the to- 
tal radiation flux (e.g. Frank, King fc Raine|2002 |: 



T(r 



where 



/ r \ -3/4 r / 7? \ 1/2' 



1/4 



(1) 



(2) 



M is the black hole mass, Af is the accretion rate and Rin 
is the inner disk radius. The emitted spectrum is calculated 
assuming that the disk emits locally as a blackbody, so that: 



Lj, = AnrgFi, = 8n cos i 



Iu{T)rdr, 



(3) 
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where i is the disk's inclination angle with respect to the 
direction of the line of sight, ro is the distance to the ob- 
server, Rout is the outer disk radius and Iv{T) is the spe- 
cific intensity. We assume that the gas is in LTE and hence 
h{T) = B^(T), where B^{T) is the Planck function. 

The MCD model is an approximation of the standard 
disk and differs from it only for the expression of the tem- 
perature profile, in which the effects of the viscous torque at 
the inner boundary are neglected. The temperature profile 
is then given by ( [Mitsuda et al.||T984| ): 



T (r) = r„ 



Rin,I 



-3/4 



where 



Tin — 



3GMM 



1/4 



(4) 



(5) 



and Rin,BB is the inner radius of the MCD. The emitted 
spectrum can again be calculated from equation Q, with 
the temperature profile given by equation Q . 

In the following, we will use the MCD model as reference 
for fitting the standard and relativistic disk spectra, because 
this model is frequently adopted for the spectral fits of X-ray 
binaries. The fitting parameters of the MCD model, the disk 
inner temperature Tin and the spectral normalization factor 
Kbb = cos i{Rin,B B / ro)'^ , depend on M and M. Assuming 
that the inner radius Rin.BB ~ bRg (Rg = GM/c? is the 
gravitational radius), it is possible to write M and M as: 



M 


67.5 / 




b \ 


M 




= O.lb^ 


MEdd 





IMpc 



1/2 



IMpc 



(KbbV^^ fT 
\ cos i J \ IkeV 



(6) 
(7) 



where Kbb = \{Ri 

„,sB/lKm)2(ro/10kpc)-2 • cos I . 
In these assumptions, if the distance ro is known and 
the inclination angle i is fixed, the MCD parameters can 
be used for estimating A4 and/or M. These expressions de- 
pend on the value of 6, the inner radius in units of Rg, and 
Tin- It is well known that these estimates of the BH mass 
and the disk accretion rate are uncertain, because the phys- 
ical state of the disk may vary and, in addition, radiative 
transfer effects induced by scattering may cause departures 
of J,y from Bv Several non-LTE accretion disk models, in 
which radiative transfer effects at the surface of the disk are 
computed, have been presented in the literature, for both 
hot disks around stellar mass BHs (e.g. [Merloni et al.|2000{ 
Davis et al.||2005 ) and cool disks around intermediate mass 
BHs (e.g. |Hui et al.|2005| . These effects are very important 
in predicting the local spectra emitted at a given disk an- 
nulus and can be accounted for using an hardening factor 
/ = Tcoi/Teff, where Tcoi and Teff are the color and effec- 
tive temperature. It is easy to show that Kbb o c l/.f* and, 



as a consequence, M cx and M /MEdd oc ( [Kubota et 
al.|1998|[200l| . So, as a first approximation, we account for 
radiative transfer effects simply inserting into equations 
([6| and ([7|. Computations of / for standard and relativistic 
disks give values in the interval 1.4-2 (seelShimura & Taka- 
hara|1995|[Zampieri et al.|2001|[Hui et al.i2005, Davis et al. 



2005b 



2.2 Relativistic accretion disk 

The structure equations and the emitted spectrum of the rel- 
ativistic disk are obtained in the same assumptions stated 
above: (1) the accretion disk is optically thick, geometri- 
cally thin and in a quasi-steady state; (2) the local emitted 
spectrum from a small annulus of the disk is a blackbody 
(/. = B^). 

In the following, we will refer to the radius of the inner- 
most stable circular orbit around a black hole with Tms- Its 
expression is given by: 



M {3 + Z2~ [(3 -Zi){3 + Zi + 2^2)]'/'} , 



(8) 



where Zi and Z2 are functions of AI and a ( [Novikov fc 
Thorne 19731, and a is the specific angular momentum of 



the black hole. 



2.2.1 Temperature profile of the relativistic disk 

A detailed derivation of the structure of a relativistic ac- 
cretion disk around a black hole is presented in [Novikov fc| 
Thorne] \1973\ and [Page fc Thorne] |l974[) to which we refer 



for all the details. Here we present only the relevant equa- 
tions for the temperature profile and the locally emitted 
flux. 



Following Novikov & Thorne (19731, the emitted radi 



ation flux can be expressed as a function of the integrated 
— r component of the stress-energy tensor (evaluated in the 
comoving frame of the observer in geodetic circular motion) 
as: 



F 



SGMM 



(9) 



where B, C and Q depend on the radial coordinate r and the 



speciflc angular momentum of the black hole a ( Novikov & 
|Thorne||1973[ ). The function Q depends also on the angular 
momentum per unit mass of a generic circular orbit and of 



the innermost stable circular orbit (Stoeger 19761. Assuming 
that the radiative transport of energy is dominant over the 
turbolent one, in LTE we can write 

F = o-T*. (10) 

Comparing equations ^ and ( 10 1 it is possible to obtain the 
surface temperature profile for a relativistic accretion disk: 



T(r) = r* 



-3/4 



1/4^1/8^1/4 



B-"^cr 



(11) 



where T« is given by equation ([2|. 



2.2.2 Relativistic effects on the propagation of radiation 



Following Bardeen et al. ( 1972 \ and Cunningham ( 1975 1 



the spectrum of a relativistic accretion disk is calculated 
by means of a trunsfer function /, that accounts for all 
relativistic effects on the emission and propagation of ra- 
diation (Doppler shift, light bending and gravitational red- 
shift). The expression for / is ( Cunningham|19"75 I: 

a ( ^ .2N1/2 2 
f (g ,re,'do)dg dre ^ (g - g' ro cos i9orf^^o, (12) 



where re is the radius of a generic emitting ring of the disk, 
t?o is the polar angle between the line of sight of the distant 
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observer and the disk's polar axis, dflo is the solid angle 
that a specific geodesies family subtends to the observer at 
a given energy and for a given emission radius, g = Eo/Ee is 
the ratio of the observed photon energy and the emitted one, 
g* , gmax and gmin are functions of g (for their expression 
see 



Cunningham|l975 1 



The radiation flux seen by a distant observer is: 



Jo cos dodTlo, 



(13) 



where /o is the observed specific intensity. If /e is the emit- 
ted specific intensity, from the exp ression of the relativistic 
invariant I / E'^ and equation ( |l3[ ) we obtain the observed 
luminosity per unit energy 



AtitqEo — 47V 



g'^ leVQ COS ^odQo- 



(14) 



Using equation (121 it is possible to rewrite tq cos-dodflo in 
terms of / {g* ,r^, j?o) dg*dr^. Inserting the resulting expres- 



sion into equation (14 1, we obtain 



2nf , 



'1/2 



dg*d (j. 



(15) 



Tabulated values of the transfer function for a Schwarzschild 
and a maximally rotating (a = 0.9981) Kerr BH were re- 



ported by Cunningham ( 1975 1. In the following we will con- 



sider these two extreme cases as reference. 




o 



log frequency (Hz) 





Disk. 






Disk. 


.stand 















1 2 3 

log B/R„ 



3 PARAMETER ESTIMATES FROM 
SPECTRAL FITS 

We implemented the structure equations and spectra of the 
accretion disk models outlined above and produced fits of 
the standard and relativistic disks against the MCD model. 
This model-model comparison has been performed by means 
of a standard minimization procedure, that provides also 
an assessment of the statistical significance of the fit. A fixed 
reference error of ~10% was assigned to the spectra. As 
mentioned above, we use the MCD model as reference for 
fitting all the other disk spectra, because this model is fre- 
quently adopted for the spectral fits of X-ray binaries within 
XSPEC. The values of the best-fitting MCD parameters are 
then compared with those of the other models. 

From now on we refer to the different accretion disk 
models with the following expressions: Disk-stand for the 



standard accretion disk ( Shakura fc Syunyaev|[l973 1 , Disk- 
BB for the MCD (iMitsudaet al.|1984||, Disk-rel or Disk-kerr 



for a relativistic disk around a Schwarzschild or a maximally 



rotating (a — 9981) Kerr black hole (Novikov & Thorne 
1973), respectively. We will also use different indeces to de- 



note the mass M and the accretion rate M of the different 
models {Mbb, Mbb for the Disk-BB; Mstand, Mstand for 
the Disk-stand; Mrei, Mrei for the Disk-rel; Mkerr, Mkerr 
for the Disk-kerr). The accretion rate is expressed in units 
of the Eddington accretion rate MEdd ~ D^ddjf? ^ where 
Ledd is the Eddington luminosity. For an accretion efficiency 
r\ = 0.1 — 0.4, the accretion luminosity reaches L^dd for 
M = MEddh = 2.5 - lOMfidd. 



Figure 1. Spectrum (top panel) and temperature profile bot- 
tom panel) of the Disk-stand model (solid-green) for M^t^nd = 
IOOM0, Mstand = O.lM^jdd, Rin = &Rg, and the correspond- 
ing best-fitting Disk-BB model (dashed-red; Mbb = 94Mq, 
Mbb = OWMEdd, Rin,BB = 13.9Rg). These plots refer to the 
fitting results reported in Table [l] second row (cosi = 1). 



3.1 Disk-stand model 

Fits of the standard disk with a MCD were performed in 
the literature in order to obtain correction factors for the 
inner disk radius, the black hole mass and the accretion 
rate (e.g. |Kubota et al.|[l998l |Makishima et al.||2000[). We 



performed similar fits of the Disk-stand with the Disk-BB 
but keeping either the BH mass or the accretion rate fixed. 
Acceptable fits with different best-fit values of Rin,BB were 
obtained (see Table [T]). For the fit with fixed M, we obtain 
Rin.BB = 13-Rg, while for the fit with fixed M, Rin^BB ~ 
14.9-Rg. In the first case, the value of Mbb returned by the 
fit is « 20% smaller than Mstand while, in the second case, 
the values of Mbb derived from the fit is « 12% smaller 
than the corresponding value of Mstand- Clearly, there are 
also fits with fixed Rin,BB in the interval 13^ Rin,BB ^14.9, 
that return intermediate values of Mbb and Mbb (see again 
Table [1]). In Figure [l] we show the results of one of these fits. 
These results are in agreement with those of |Kubota et| 



al 



1998[ ), who showed that the inner radius of the best- 
fitting MCD model does not correspond to the true inner 
radius of the standard disk, but that the two radii are related 
by the expression (see also Zimmerman et al.|2005 I 
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Figure 2. Same as Figure [T] for tlie Disk-rel model (solid- 
orange; Mrel = WOMq, Mrel = O-lMEdd, R,-n,rel = GRg). The 
parameters of tlie best-fitting Disk-BB model (dashed-red) are 



MbB = 87Mq, MbB = O.OYMEdd, Rrr,,BB 

ble[2] fifth row). For both models cosi = 1. 



2l.&Rg (see Ta- 



Figure 3. Same as Figure [T] for the Disk-kerr model {solid- 
blue; Mi^^^j. = 1OOM0, Mf^g^^j. = O.lMBdd)- The parameters of 
the best-fitting Disk-BB model (dashed-red) are Mbb = 88Mq, 

Mbb = 0.095MEdd, Rin,BB " 
For both models cosi = 1. 



8.9Rg (see Table [s] second row). 



Rin,BB — 2.6i?i„ 



(16) 



In fact, in order for the spectral peaks of the Disk-BB 
and Disk-stand be superimposed (see Figure the max- 
imum temperature of the standard disk (0.488r»; see equa- 
tions [l] and [2]) must be equal to the maximum temper- 
ature of the Disk-BB (Ti„; see equation [5]). From this, 
equation ( |16[ ) easily follows. Therefore, assuming that the 
Disk-stand terminates at the ISCO (innermost stable circu- 
lar orbit), Rin = 6Rg and the inner radius of the best-fitting 
Disk-BB must be located at ~ 15.6Rg, in subtantial agree- 
ment with the results of the model-model fits. 

We note that all the results presented here hold for any 
inclination angle, as the Disk-BB and Disk-stand models 
have the same dependence on i. 



3.2 Disk-rel model 

Similar fits of the Disk-rel model with the Disk-BB were also 
performed. For sufficiently small inclination angles, also in 
this case the emitted spectrum can be well reproduced by 
the MCD. The fits were performed keeping either M or AI 
fixed. Results are shown in Table (2] The two fits return 
different values of the inner radius, i?m,ss ~ 19.2Rg for 



fixed M and Rin,BB 
the value of Mbb returned by the fit is ~ 40% smaller than 
Mrei wheras, in the second case, the value of Mbb derived 
from the fit is « 25% smaller than the corresponding value 
of Mrd- From Table [2] we can see that good fits can be 
obtained also fixing Rin,BB in the interval 19.2-25i?g and 
leaving Mb b and Mb b free to vary. In Figure 1] we plot the 
results of one of these fits {Rin,BB ~ 21.6Rg). 

If we increase the inclination angle, the relativistic ef- 
fects (especially the gravitational focusing) tend to increase 
the contribution of radiation coming from the inner parts of 
the disk, making the spectrum of the Disk-rel intrinsically 
harder than that of the Disk-BB. In fact, for high inclina- 
tion angles the Disk-BB is only in rough agreement (within 
a ~ 10% fractional error) with the Disk-rel spectrum, and 
no satisfactory fit is obtained (Xred ^ 1-3). Similar results 
for cosi < 1 were obtained by Ebisawa et al. (19911. For 
illustrative purposes alone, in Table [2| we report also the 
values of the fits for cosi = 0.5. 



3.3 Disk-kerr model 

The last comparison concerns the Disk-kerr/ Disk- BB mod- 
els. If the inclination angle is close to i = (face-on disk), 
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Table 1. Values of the characteristic parameters of the Disk-stand model and those derived by the spectral 
fits with the Disk-BB model for cosi = 1. The inner radius of the Disk-stand model is &Rg. In the first row 
Mbb is set equal to M^tandt in the second row Rin,BB is fixed at an intermediate value between 13 and 
14.9 Rg; in the last row Mbb is set equal to M^i^nd- The results of the spectral fit at fixed radius are shown 
in Figure^ 



Mstand 


Mgtand 


Mbb 


Mbb 


Rin.BB 






(Mq) 


(MEdd) 


(MEdd) 


(Rg) 




100.0 


0.1 


100.0 


O.OSitO.Ol 


13. Oil. 5 


0.1 


100.0 


0.1 


94.0±10.0 


0.09±0.03 


13.9 


0.2 


100.0 


0.1 


88.0±32.0 


0.1 


14.9±1.3 


0.2 



Table 2. Values of the characteristic parameters of the Disk-rel model and those derived by the spectral fits 
with the Disk-BB model for cosi = 1 and cosi = 0.5. The inner radius of the Disk-rel is Ri„ = 6Rg. In 
the first row Mbb is set equal to MreV, in the second row Rin,BB is fixed at an intermediate value between 
19.2 and 25 Rg; in the last row Mbb is set equal to M^-^i (for each value of cosi). The results of the spectral 
fit at fixed radius are shown in Figure [2] 
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0.1 
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0.07±0.02 
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100.0 
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100.0 
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100.0 
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0.1 
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1.8 



the spectral fits of the Disk-kerr with the Disk-BB at con- 
stant M or M are statistically acceptable and return dif- 
ferent values of the inner radius, Rin,BB = 7.9Rg and 
Rin,BB = 9.2Rg, respectively (see Table |3|. The values of 
Mbb and Mbb are larger than the corresponding values of 
the Disk-rel model (as high as ~ 20% for M and ~ 25% for 
M). Also in this case, there are other fits with fixed Rin,BB 
in the range 7.9-9.2i?g for which the Disk-BB and Disk-kerr 
spectra are in good agreement (see again Table [3|. In Fig- 
ure [3] we show the results of one of these fits obtained for 

Rin,BB ~ 8.9i?g. 

We have also tried fits varying the inclination angle i. 
For a maximally rotating Kerr black hole, the relativistic 
effects on the emitted spectrum are significantly stronger 
than for a Schwarzschild black hole. The spectrum of the 
Disk-kerr model becomes intrinsecally different and much 
harder than the Disk-BB one as i increases, so that it is 
no longer possible to obtain a satisfactory fit (x^ > 2). For 
illustrative purposes alone, in Table [S] we report the values 
of the fits for cosi = 0.5. 



3.4 Comparison with fitting results from XSPEC 

In order to test our approach, we simulated disk spectra for 
both a Schwarzschild BH and a maximally rotating Kerr 



BH using the kerrBB model in XSPEC (v. 12.3.0). This 
model implements the relativistic structure equations and 
the relativistic effects on the emission and propagation of 
radiation ( Li et al.|2005 1 . We fitted the kerrBB spectra with 
a diskBB model (the XSPEC implementation of the MCD) 
and used the parameters of the fit to estimate the mass and 
accretion rate using equations ^ and ([7|. 

Spectra were simulated using the response matrices of 
the EPIC-pn camera onboard XMM-Newton. For the sake of 
comparison, in the kerrBB model we adopted a zero torque 
inner boundary condition and switched off the effects of self- 
irradiation and limb-darkening. The disk was assumed to be 
face-on {i = 0) and the color correction factor / was set 
equal to 1. A fixed distance of 5 Mpc was assumed. For the 
mass of the BH and the accretion rate we chose lOOAf© and 
0.485M_Bdd, respectively, while the BH angular momentum 
a was set equal to (Schwarzschild) or 0.9981 (Kerr). 

For a kerrBB spectrum with a = 0, the best fit- 
ting MCD parameters are: Tin = 0.122 ± 0.002 keV and 
A' = 39 ± 6 (2(T errors). From these parameters, assuming 
b = 19.2±3.5 (see Table[2|, we obtain Mbb = (110±29)Mo 
(eq. [g]), in good agreement with the assumed value of the 
BH mass. Similarly, assuming fe = 25 ± 6.7, it is Mbb ~ 
(0.4 ± 0.28)MEdd (eq. [t]), which is also in agreement with 
the assumed value for the accretion rate. For a kerrBB spec- 
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Table 3. Values of the characteristic parameters of the Disk-kerr model and those derived by the spectral 
fits with the Disk-BB model for cosi = 1 and cosi = 0.5. The inner radius of the Disk-kerr model has been 
calculated using equation ([8|l (for a = 0.9981, Rin,kerr — l.SG-Rg). In the first row Mbb is set equal to 



Mr, 



J.; in the second row Rin.BB is fixed at an intermediate value between 7.9 and 9.2 i?„; in the last row 



Mbb is set equal to Mfcerr 
in Figure [3] 



(for each value of cosi). The results of the spectral fit at fixed radius are shown 



^Kerr 




Mbb 


Mbb 


Rin.BB 


^red 


(Me) 


(Mq) 


{MEdd) 


{MEdd) 


{R9) 




cosi = l 100.0 


0.1 


100.0 


0.075±0.015 


7.9±0.9 


0.2 


100.0 


0.1 


88.0±18.0 


0.095±0.02 


8.9 


0.2 


100.0 


0.1 


85.0±20.0 


0.1 


9.2±1.3 


0.2 


cosi = 0.5 100.0 


0.1 


100.0 


0.04±0.015 


2.4±1.7 


2.3 


100.0 


0.1 


60.0±5.0 


0.1 


4.0±0.4 


2.2 



Table 4. Values of the BH mass M and accretion rate M for NGC 1313 X-1 and NGC 1313 X-2, obtained 
using the parameters of the MCD that best-fits their soft components {Kbb = 28^g, Ti„ = 0.22 keV 
for NGC 1313 X-1 ; Kbb = 6.66+^®,,, Ti„ = 0.25 keV for NGC 1313 X-2) and different values of b that 



correspond to the different accretion disk models (see Tables [2] [3] and section 4; [Miller et al.|2003| and |Miller| 
|et al. 11200411 adopted b = 9. 5). The assumed distance of NGC 1313 is 3.7 Mpc ||Tully||1988|l and the color 
correction factor is / = 1.7 | |Shimura &: Takahara|1995| [Zampieri et al.|2001| [Hui et a.l'|2005| . 





NGC 1313 X-1 




NGC 1313 X-2 




Rin,BB {Rg) 


M (Mq) 


M {MEdd) 


M (Mq) 


M {MEdd) 


6 


630^50 


n ^-1-0.01 


3-,0+260 
'^I'-'-lOO 


4+0-3 
'-'■^-0.1 


7.9 
9.2 






24Q+190 


9+0-8 


9.5 








1-^-0.3 


13 
14.9 


290^30 


og + 0.2 


145l-« 


2 4+2-0 


19.2 
25 


200^20 


q+0.6 


lool- 


°-'-2.0 



trum with a = 0.9981, the best fitting MCD parameters 
are: Ti„ = 0.248 ± 0.003 keV and = 6 ± 0.4 (2cr er- 
rors). From equations ([6| and ([7| and the values of b re- 
ported in Table [3] we obtain Mbb = (105 ± \^)Mq and 
Mbb = (0.4 ± 0.14)M_Bdd, again in satisfactory agreement 
with the parameters assumed in simulating the spectrum. 



4 ESTIMATES OF M AND M FROM THE 
X-RAY SOFT COMPONENT OF ULXS 

A certain number of ULXs show two component X-ray spec- 
tra that can be adequately modelled with a soft thermal 
component plus a power-law. If the soft component is fitted 
with a MCD, the fitting parameters (in particular the nor- 
malization Kbb) return very large values for the BH mass, 
often in excess of several hundreds Mq (e.. 



Miller et al. (20031 



Rg) and the color correction factor /. 
and [Miller et al.| ( |2004"a| ) adopted b = 9.5 and / = 1.7. 

In the previous Section, we performed fits of the stan- 
dard and relativistic disk with a MCD keeping the BH mass 
fixed and determined the (range of) values of b for which 
the spectra are in statistical agreement. These values are re- 
ported in Tables [l]j2]and [3] and differ from what assumed 

p004al 



by Miller et al. 



standard disk. 



m 

p003 



2003| ) and [Miller et al 
there appears 



to 



For the 

be more agreement with 



what reported by [Kubota et al.[ ( [1998[ ) and [Makishima et[ 



al 



(20001 



Miller et al. 

[20031 [2004a[ ). These masses are calculated from equation (|6|) 
assuming specific values for the inner radius 6 (in units of and NGC 4559 X-7 obtained by Miller et al. (2003 2004a I, 



We are then in the position to revise previous es- 
timates of the BH mass (and accretion rate) of ULXs based 
on MCD spectral fits of the soft component, adopting the 
appropriate value of b for both a standard and a relativistic 
disk. 

For illustrative purposes, we use the results of the MCD 
spectral fits of NGC 1313 X-1, NGC 1313 X-2, M 81 X-9 
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Table 5. Values of the BH mass M and accretion rate M for M81 X-9 and NGC 4559 X-7, obtained using 
the parameters of the MOD that best-fits their soft components [Kbb = 20tio' "^'i = 0-26 keV for M8f X-9 
(obs. 1); Kbb = 60l™, Ti„ = 0.21 keV for M81 X-9 (obs. 2); Kbb = ISS^^o?. Ti„ = 0.148 keV for NGC 
4559 X-7) and different values of b that correspond to the different accretion disk models (see Tables [5] |3] 
and section 4; |Miller et al.|2003|and|Miller et al.|2004a| adopted b = 9.5). The assumed distance is 9.69 Mpc 
for NGC 4559 X-7 ( [Sanders et ar]|2003|l and 3.4 Mpc for M81 X-9 ||Georgiev et al.||199l"] |Hill||1993| |. The 
color correction factor is / = 1.7 jShimura &: Takahara|1995| [Zampieri et al.|2001| |Hui et al.|2005[ l. 





M 81 X-9 (obs. 1) 




M 81 X-9 (obs. 2) 




NGC 4559 X-7 
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M (Mq) 
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[Zampieri et al.|(2004| ) and |Cropper et al.[ ( |2004[ ). From equa- 
tions ^ and ([7 1 of Sectioii ]"2.1[ we calculate the values of M 
and M for different values of 6 and assuming / — 1.7, and 
report them in Tables |4] and |5] The values that correspond 
to 6 = 9.5 are those reported by Miller et al. (20031 and 
iller et ah] ( |2004a[ ) , whereas the values for b = 7.9, 13, 19.2 
correspond to those expected for a standard disk around a 
maximally rotating Kerr, newtonian, or Schwarzschild BH, 
respectively. As can be seen from the Tables, the values of 
the mass are systematically lower (up to a factor 3) than 
those of a MOD with R^u.bb at the ISCO (& = 6). This is 
clearly a consequence of the larger fitting radius required by 
the newtonian and relativistic disks. Comparing the masses 
of the different "physical" models, it is possible to see that 
the largest value is obtained for a disk around a maximally 
rotating Kerr BH [Disk-kerr; see Tables [4] and [5| . This is 
consistent with the fact that the (Boyer-Lindquist) radial 
coordinate of the ISCO for a maximally rotating, Kerr BH 
is ~ 1/6 that of a Schwarzschild BH of the same mass (cmp. 



Hui & Krolik 20081. Indeed, for a fixed value of Rin,BB 



(in cm) as returned by the MCD fit of an observed spec- 
trum, a lower value of 6 necessarily implies a larger mass as 
b — Rin,BB / Rg oc 1/AI. We uccd to fill the same physical 
radius Rin,BB with a bigger BH in order to reduce the value 
of b. Ultimately, the smaller value of 6 returned by the fit 
of the Disk-kerr model is a consequence of the fact that the 
Kerr spectrum is significantly harder than those of the other 
"physical" disk models computed with the same value of the 
parameters (see Figure |4|. 

It is also interesting to note that, apart from NGC 4559 
X-7, the BH masses for a Schwarzschild disk are no longer 
significantly in excess of ~ IOOMq, with a distribution clus- 
tering in the interval 100-200 Mq (see again Tables |4] and [5|. 
This is true also for M 81 X-9, for which the second observa- 
tion returns somewhat larger values of the mass. However, as 



the limits inferred from the two observations of this source 
must be consistent within the errors, the mass of the BH 
must be in the interval 150 — 220Mq (smaller/larger values 
would not be consistent with the second/first observation; 
see the last row in Table [5|. All the mass estimates obtained 
for the Disk-rel model are then significantly lower than those 
reported by [Miller et"ar]|2003j and Miller et al.[ ( [2004a| . The 
only exception is NGC 4559 X-7 for which, however, there 
are hints of timing features in the power density spectrum 
and it has been proposed that it may contain a BH of sev- 
eral hundreds solar masses (e.g. Cropper et al. 20041. The 



inferred BH masses can reach the values estimated by Miller[ 
[et al.[ ( [2004a[ ) only for a disk around a maximally rotating, 
Kerr BH. 



5 CONCLUSIONS 

The spectrum of a standard disk can be well reproduced by 
a MCD model. However, as already noted by [Kubota et al.] 
(19981 [2001| [2005|, if one tries to exploit the MCD model 



to derive the values of the accretion rate and black hole 
mass, with the usual assumption that, in the soft state, the 
disk terminates at the ISCO {Rin,BB = 6J?g), the inferred 
values of AI and AI are largely overestimated. On the other 
hand, the spectrum of the standard disk can be fitted by 
a MCD model with the same values of M or M if we take 
Rin,BB = bRg with 6 ~ 13 or & ~ 15, respectively. This is in 
substantial agreement with the approximate value (6 = 15.6) 
that can be derived using analytic arguments. 

Also in the case of a relativistic disk around a 
Schwarzschild black hole, for sufficiently small inclination 
angles the emitted spectrum can be well reproduced by a 



MCD model (see also Ebisawa et al. 1991 and Kubota et 
[al. 12005] ). For small angles, in fact, one can try to guess the 
disk and black hole parameters using the best-fitting MCD 
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log frequency (Hz) 




space in which the MCD is in agreement with "more phys- 
ical" standard accretion disk models and provide "recipes" 
for adjusting the estimates of the disk inner radius, the black 
hole mass and the accretion rate. In fact, despite the un- 
certainties related to radiative transfer and to the physical 
state of the inner disk, the parameters of the MCD fit are 
often used for inferring the mass and accretion rate in X-ray 
binaries, in particular in ULXs. The results of the compar- 
ison that we performed can be used to revise estimates of 
AI and M obtained in this way. This has also been tested 
by applying our "recipe procedure" to the results from of a 
diskBB fit to simulated kerrBB spectra performed directly 
in XSPEC. 

We considered the case of a few ULXs (NGC 1313 X- 
1, NGC 1313 X-2, M81 X-9 and NGC 4559 X-7) for which 
MCD spectral fits of their X-ray soft spectral components 
have been published and/or values of the BH mass estimated 



(Cropper et al. 2004 Miller et al. 2003 2004a Zampieri et 



al. 20041. From the parameters of the fit with the Disk- 



los R/R„ 



BB {Tin and Kbb) we found that, assuming that the inner 
disk boundary of the MCD is at or close to QRg, the values 
of M and M can be severely overestimated. Adopting the 
appropriate value of the inner radius b (in units of Rg) for 
the Disk-rel model, we obtain that the BH masses are in 
the range ^ 100 - 2OOM0 for NGC 1313 X-1, NGC 1313 
X-2, M81 X-9, and 700 - 2OOOM0 for NGC 4559 X-7 (see 
Tables [4] and |5|. The relativistic disk models have BH masses 
systematically lower (up to a factor 3) than those of a MCD 
with Rin,BB at the ISCO. 



Figure 4. Emitted spectra (top panel) and temperature profiles 
{bottom panel) of the Disk-BB, Disk-stand, Disk-rel and Disk-kerr 
models, obtained for the same characteristic parameters: M = 



As already mentioned, sophisticated relativistic accre- 



IOOMq, M = O.lMEdd, Rin,BB = Rin 



R, 



Rin,kerr has been obtained using equation For all models 
cosi = 1. 



model but, as for the standard disk, a suitable assumption 
must be made on the inner disk radius. We found that there 
exist precise values of Rin,BB {— 19.2Rg and ~ 25Rg) for 
which the values of M and M inferred from the MCD are 
the same as those provided by the relativistic disk. For large 
angles, owing to genuine relativistic effects, no satisfactory 
MCD spectral fit of the relativitisc disk spectrum can be 



obtained (e.g. Ebisawa et al.||1991 1 



The spectrum of an accretion disk around a fast rotat- 
ing Kerr black hole, is in general more strongly affected by 
relativistic effects compared with those of a disk around a 
Schwarzschild blak hole. As a consequence, it can be well re- 
produced by a MCD model, only in the case in which the disk 
is essentially face-on. In this assumption, there exist values 
of Rin,BB for which the MCD and Kerr disks return similar 
values of M and M {Rin,BB — 7.9Rg and Rin,BB — 
respectively). Also in this case, for large inclination angles, 
the spectrum of the Disk-kerr model becomes intrinsecally 
different and much harder than the Disk-BB one, so that 
it is no longer possible to obtain a satisfactory fit with the 
MCD model. 

The model-model comparison approach adopted in this 
Paper was prompted by the idea to determine the parameter 



tion disk models (e.g. Hanawa||1989 Li et al. 2005 Hui et 



al.|2005 \ that implement the relativistic structure equations 



and the relativistic effects on the emission and propagation 
of radiation, are available in XSPEC and can clearly be used 
to perform accurate fits of observed spectra, in particular for 
large inclination angles when the relativistic effects produce 
significant differences. An analysis of this type has recently 



been performed by Winter et al. (20071 and Hui & Krolik 



( |2008[ ). It is interesting to note that the BH masses for a 
sample of disk-dominated ULXs obtained by |Hui fc Krolik| 
(200^ are below ~ IOOMq. Even if their best fits require 
large values of the BH spin, BHs of several hundreds solar 
masses are not needed, in substantial agreement with our 
findings. 

Clearly, it is possible that at least some ULXs for which 
curvature above ~ 2-3 keV is present in the X-ray spectrum 
are in a different accretion regime, in which the disk is in a 
diverse physical state (accreting at or above the Eddington 
rate; see e.g. |Roberts|[2007[ |Soria||2007l ). In this case, the 
assumptions of a standard accretion disk break down and 



different spectral models should be adopted (e.g. Stobbart 
|et al.|2006||Mizuno et al.|2007| . However, in the assumption 
that they are in a standard accretion regime (below Edding- 
ton), we found that, in three cases (NGC 1313 X-1, X-2 and 
M 81 X-9), the black hole masses inferred for a standard 
disk around a Schwarzschild black hole are in the interval 
~ 100 — 2OOM0. Only if the black hole is maximally rotat- 
ing are the masses comparable to the much larger values 
previously derived in the literature. 
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